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ABSTRACT 



Esmailzadeh et al 



(I1991I ) to calculate 



We use the semi-analytic method of 
the abundances of Helium and Deuterium produced during Big Bang nucleosyn- 
thesis assuming the fine structure constant and the Higgs vacuum expectation 
value may vary in time. We analyze the dependence on the fundamental con- 
stants of the nucleon mass, nuclear binding energies and cross sections involved 
in the calculation of the abundances. Unlike previous works, we do not assume 
the chiral limit of QCD. Rather, we take into account the quark masses and 
consider the one-pion exchange potential, within perturbation theory, for the 
proton-neutron scattering. However, we do not consider the time variation of the 
strong interactions scale but attribute the changes in the quark masses to the 
temporal variation of the Higgs vacuum expectation value. Using the observa- 
tional data of the helium and deuterium, we put constraints on the variation of 
the fundamental constants between the time of nucleosynthesis and the present 
time. 



1. Introduction 

Big Bang Nucleosynthesis (BBN) offers the deepest reliable probe of the early universe. 
Predictions of the abundances of the light elements D,^He, "^He and ^Li synthesized at the 
end of the 'first three minutes' are in good overall agreement with the primordial abundances 
inferred from observational data, which validates the Standard Big Bang Nucleosynthesis 
(SBBN). BBN also provides powerful constraints on possible deviations from th e stand ard 



cosmology and on new theories on physics beyond the Standard Model (SM) (ISarkar 



19961 ). Among these theories, there are those in which some of the dimen sionless ratios 



of fundamental constants do vary in time like string-derived field theories (IWu and Wang 
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1986 



2002a 



1983 



Maedalll988l : iBarr and Mohapatra Il988l: 



b|), related brane- world theories (lYoum 



D amour and Polvakov 



2001a 



2003 ) , and (related or not) K a. 



Gleiser and Tavloi 



astronomical data ( jWebb et al. 



1985 



uza-Klein theories (Kaluza 



b 



Pa 



Overduin and Wesson 



1999 



2001 



Murphy et al 



ma et al 



1921 



199 



2003 



Klein 



Damour et al. 



Brax et al.l 



1926 : 



Weinberg 



1997 ) . On the othe r hand, recent 



2001a b 



20031 ) suggest a 



possible variation of the fine structure constant a = e^/hc at the 10 ^ level over a 



time-scale of 10 billion years. However, other recent 



( Martinez Fiorenzano et al. 



2004 



Grupe et al. 



200^ 



2003 



Quast et al. 



Chand et al. 



measurements of molecular hydrogen (jivanchik et al. 



2004 



independent analys i s of similar dat a 



Bahcall et al 



2004 



Srianand et al.l 



20061) found no var i ation . On the other hand 



2002, 



2003 



20051 ) reported a variation 



of the proton to electron mass /i = — . This fact motivated more gene ral discussions o 



possible variations of other constants. 



Calmet and Fritzschl (120021 ) and 



Langacker et al. 



( 2002) have studied the implication of gauge unification for the time variation of a while 



Olive et al. 



(IBekenstein 



2OO2I) explored a super-symmetric version of the dynamical Bekenstein model 



19821 ) in order to produce a large change in a in the redshift range z = 0.5 — 3.5 



and still be consistent with the constraints on Aa/a from the results of high precision 
limits on the violation of equivalence principle by a fifth force. 

On the other hand, there are many non-SBBN models which introduce new free 
parameters in addition to the baryon density parameter, or equivalently the baryon 
asymmetry rjB = = 2.74 x lO^^Qsh"^- Most known of these models are those which 

assume either a non-standard contribution to the total density, or a lepton asymmetry. The 
first possibility affects the expansion rate of the universe S = ^ = and can be restated 
in terms of 'equivalent' number of extra neutrinos AN^ = — 3. Simple analytic fits to 
BBN and the cosmic mi crowave background ra diation (CMBR) dat a provide the following 



bound: 0.85 < 5 < 1.15 ( IBarger et al. 



2003a 



b 



Steigman 



2005, 



20061 ) . As regards the lepton 



asymmetry, observational data do not imply that is should be connected to the 'tiny' baryon 
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asymmetry rjE, and it could be large enough to perturbe SBBN predictions . Moreover, a 
small asymmetry between electron type neutrinos and antineutrinos can have a significant 
impact on BBN since the z/g affect the interconversion of neutrons to protons changing 
the equilibrium neutron-to-proton ratio from {n/p)^^ = to {n/p)eq = {n / p)^^^e~^''' . In 

consequence, the "^He abundance changes. In contrast the D abundance is insensitive to 
7^ 0. Consis tent with the BBN and CMBR data , value s of in the range —0.1 < < 0.3 



are permitted (iBarger et al 



2003a 



Steigman 



2005 



20061 ). In our analysis, however, we shall 



not consider these non-SBBN scenarios, but attribute any non-SBBN issue to time- variation 
of fundamental constants. 

The de nsity of baryonic matter Q r/i^ can be estimated using the WMAP data from 



the CMBR f lSpergel et al. 



2003 



20061 ). From the observed WMAP baryon density, the 
predicted abundances are highly consistent with the observed D but not with ^He and ''Li. 
However, any change in the value of the fundamental constants would work its way into 
the value of the abundances of the various light elements and the question we address is 
whether or not existing observations of the primordial abundances suggest any change in 
the values of the fundamental constants at the time of BBN. 

BBN is sensitive to a number of fundamental dimensionless parameters including 
the fine structure constant a , I^qcdI Mpiank and mfl j Kqcd where m'^ is the quark mass 
and IS.QCD is the strong scale determined by the position of the pole in the perturbative 
QCD running coupling constant. Several authors have studied the dependencies of 
the BBN parameters on the fundamental constants. The depend ence of the prim o rdial 



abundances on t he fine structure co n stant h as been evaluated W 



and improved by 



NoUett and LopezI (120021 ). lYoo and Scherren (120031 ) analyze the effects 



Bergstrom et al. 



(119991 ) 



of the time variation of the higgs vacuum expectation value < f > on BBN and th e 



cosmic microwave background radiation (CMBR). On the other hand. 



Miiller et al 
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( 120041 ) calculate the abundances as function of the Planck Mass Mp, a, < f >, electron 
mass rrie, nucleon decay time, deuterium binding energy {eo) and neutron-proton mass 
difference (Am = m„ — rrip). Moreover, they study the depend ence of the last three 



quant ities as functions of the fundamental couplings and masses. 



Kneller and McLaughlin 



( 20031) study the dependence of the primordial abundances on en and Am . These papers 



(lYoo and Scherreri 12003 



Miiller et al 



of chiral perturbation theory (IBeane and Savage 



2004 



Kne 



ler and McLaughlin 



20031 ) use the results 



20031 ) to adress the dependence of the 



deuterium binding energy with the higgs vacuum expectations value and Aqqd- However, 
the dependence of the deuterium binding energy on < f > was estimated from an 
approximated linear dependence of e/j on the pion mass m^r, while the exact limits on the 
relati ve change of < v > would d epen d on the details of such depend ence. On the other 



hand, 



Campbell and Olivd (Il995l ) and 



Ichikawa and Kawasakil (120021 ) study the effects 



of variation of fundamental constants on BBN in the context of a dilaton superstring 
model. Finally, limits on cosmological variations of a, Aqcd and m'^ from optical quasar 



absorption spectra, labo r atory atomic clocks anc 



Flambaum and Shuryaki (120021 ) 



Flambaum et al 



Tom BBN have been established by 



(120041 ). For computing the deuterium 



binding energy (e/) ~ 2.225MeV) they apply quantum mechanics perturbation theory. This 
factor is very significant in influencing the reaction rate of p + n—^d + '-f which is the first 
and most crucial step in BBN. 



The BBN abund ances can be computed 



19921 ) and analytica . 



work (iLandau et al. 



Esmailzadeh et al. 



u sing numerical (Wagoner 



1991 



Mukhanov 



1973 



Kawano 



20031 ) methods. In a previous 



20061 ). we used the semi-analytic method of lEsmailzadeh et al.l (119911 ) 



to calculate the abundances of the light elements produced during BBN assuming that the 
gauge coupling constants may vary in time. We considered the chiral limit of QCD when 
analyzing the nucleon masses, binding energies and the cross sections. Deviations between 
standard cosmology calculations and observational data could be interpreted as resulting 
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from variations in Gp the Fermi constant, a the fine structure constant and Aqcd the 
strong interactions scale. The semi-analytical method allows us to obtain semi-analytic 
dependencies of the primordial abundances on the fundamental constants, which otherwise 
must be computed using numerical codes. 



On the other hand, in the standard model, a variation of the Fermi constant implie s a 



variation of the vacuum expectation value (vev) of the Higgs field ( iDixit and Sher 



19881): 



^ ~ V2M^ ~ V2<v>^ 
Here Mw is the mass of the W-boson, and < v >~ 250GeV is the vev of the Higgs field. 

Within the QCD chiral limit, the quark masses, which are also proportional to the Higgs 

vev through the relation m'^ = Yyukawa < v >, are neglected. Therefore, the logical and 

consistent step to follow is to go beyond the chiral limit and take the variation of < f > 

as affecting the Fermi constant as well as the quark masses. We analyze the nucleon 

masses, the nuclear binding energies and the cross sections dependence within quantum 

mechanics perturbation theory. The objective of this paper is to study such variations as 

model independent as possible. Therefore, we consider the one-pion exchange potential 

as the perturbation on the p-n scattering responsible for the formation of the deuterium. 

This perturbation potential varies in time if the pion mass changes in time which leads 

to a time variation of the deuterium binding energy. The pion mass also depends on the 

Higgs vev through the Gell-Mann-Oakes-Renner relation: = m''i^||^ ~ tti'^Aqcd, 

where < gg > is the quark condensate and is the pion decay constant. In order to 

determine the dependence of the deuterium binding energy on the fundamental constants, 

we use the square well model to approximate the attractive strong interaction potential 

of the deuterium and fit current scattering data to get estimates for the depth and width 

of the well. On the other hand, we will not discuss the effect of Aqcd variation on the 
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QCD-determined quantities such as the quark condensate or the width and depth of the 
square well. The reason for this is that we lack a complete theory for these quantities, 
and, especially, because of the absence of p-n scattering data in the far past. For the same 
reason, we will not consider changes in the Yukawa couplings either. On the other hand, 
the effect of of A qcd variation on the abundances of the light elements was analyzed in a 



previous work (iLandau et al. 



20061 ). Thus, we will limit ourselves in this paper to study 



the dependence on a and < f > for the physical quantities, such as binding energies, 
nucleon masses and cross sections involved in the BBN calculations. Our treatment of the 



deuter ium binding energy is similar to the one performed by of Flambaum and Shuryak 



( I2OO2I ). However, there are some technical differences in the wave function normalization 
which we describe in section |H Furthermore, we go one step further in calculating the 
effects of such variations in the fundamental constants on the primordial abundances of D 
and '^He. This is one of the advantages of using the semi-analytical method and it also 
allows us to compare with observational data in order to put bounds on the variation of a 
and < V >. 

On the other hand, the concordance between the WMAP estimates and SBBN has 
been investigated 



2004 



Cyburt 



) y many authors (^Cvburt et al. 



2004 



Coc et al. 



2004a 



2003 



Romano et al. 



2003 



Cuoco et al. 



b|). From the WMAP baryon density , the predicted 
abundances are highly consistent with the observed D but not with ^He and ''Li. They 
are produced more than observed. Such discrepancy is usually ascribed to non reported 
systematic errors in the observations of ^He and ''Li. Indeed, more realistic determ inations 



of the ^He uncertainty irnplies a baryon density in li ne with the WMAP estimate fi Cyburt 



2004 



Olive and Skillman 



20041 ). On the other hand. 



Richard et al. 



(I2OO5I ) have pointed out 



that a better understanding of turbulent transport in the radiative zones of the stars is 
needed for a better determination of the ^Li abundance. In our previous work, we obtained 
results consistent with variation of fundamental constants when considering all data. 
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However, discarding the ^Li data we obtained no variation. Therefore, we suspect that 
the possible non reported systematic uncertainties are "hidden" within a setup involving 
variation of the fundamental constants. Thus, until better estimations of the systematic 
errors of ^Li are reported, we will only consider the D and ^He data. 

Even though the WMAP estimate of the baryon density is the most accurate one, it 



is sti ll affected by degeneracies with other cosmological parameters (jSpergel et al 



2003 



20061 ). On the other hand, this quantity can be also determined combining data from 
galaxy surveys (SPS S, 2dF) and x-ray satellites (Chandra, XMM- Newton, ROSAT, ASCA) 



( iLandau et a 



estimate and 



20061). In th i s wor k, we consider a weighted mean between the WMAP 



Landau et al. 



(120061 ) estimate for f2f,/i^, and, furthermore, we shall compute 
the dependence of binding energies, cross sections and abundances on this parameter. 
Finally, we shall use observational data from D and ^He to estimate the variations in time 
of a and < v > and a possible deviation of from its considered value. We also compare 
our results with other non-SBBN models, where a non-standard expansion rate and an 
electron-neutrino asymmetry were considered. Finally, we would like to emphasize that the 
approach in this work is phenomeno logical and the results we get are model independent. 

The paper is organized as follows. In section 2 we present the notations used and 
summarize the steps which one follows in the semi-analytic approach to calculate the 
abundances. In section 3 we calculate the dependence of the abundances on a, < f >, 
Ofo/i^ and the deuterium binding energy e^. In section 4, we express the dependence of 
the deuterium binding energy on the Higgs vev within the square well model. Results of 
comparing theoretical prediction with observational values are presented in section 5, where 
we also compare with other non-standard BBN models results. Conclusions are presented 
in section 6. 
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Preliminaries 



The method of lEsmailzadeh et aLl (jl99ll ) consists of calculating the different 
abundances between fixed points or stages. One solves the equations for the light elements 
only for one element in each stage. For the other elements (say, the i^^), it is necessary to 
solve the quasi static equilibrium (QSE) equation {Yi ^ 0), where Yi is the abundance of the 
jth gig]2ient relatively to baryons, considering only the most important rates of production 
and destruction. On the other hand, we should also calculate the final temperature of 
each stage. We show in Table [1] the different stages and their corresponding equations, to 
which is added als o the conservation of t he ne utron number (further details are given in 



Esmailzadeh et al. 



(|l99lh 



Landau et al. 



mod }). 



Table 1: Stages and equations, n refers to neutron, p to proton, d to deuterium, T to tritium, 
3 to •^He and a to ^He 



Stage 


Equations 


Final temperature 


Until the weak interaction freeze-out 






Until the production of ^He becomes efficient 


Y^ = -2Y^-Y^[n] 
Y, = Y3 = Yt = 


2Y^ ~ Y^[n] 


Until the production of deuterium dominates 
rate of change of neutrons 


Y — —2Y 
Yd = Ys = YT = 


Yn = Ya 


Deuterium final abundance 


Yd = -2Y^ 

Yn = Ys = YT = 





Since we are considering changes in the Higgs vev (< f >) and the fine structure 
constant (a), we need to find expressions for the nucleon masses and binding energies in 
terms of these quantities. For the P-N mass splitting we have 
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Q = Am = rrin — rrip = A"m + 
where A"m is the contribution of the electromagnetic energy, and thus 



(2) 



^. On the 



other hand, is due to p-w mesons mixing and known to be proportional to 



( IChristiansen et al. 



1991 



Therefore, 



5mi 
ml 



S<v> 



6Q 
Q 



<v> 

6a 6 < V > 
-0.587— + 1.587 

a < V > 



. Thus, we get: 



(3) 



We need also to know the dependence of the variation of the nuclear mass of an 
element iii terms of the changes in < f > and a. This can be estimated using 
M{X) = ZrUp + Nnin — ex where ex is the binding energy for the element X, and we have: 



Se 



X 



ec Sa 
ex a 



(4) 



where ec 



Z 



is the electromagnetic contribution. The radius of the nucleus 



{R ~ 1.2Aifm) is considered as a strong interaction effect and, thus, taken to be constant 
in our analysis. The change in the neutron decay rate in terms of the changes in a and 



< V > can be expressed as follows (jlchikawa and Kawasaki 



2002 



Landau et al. 



20061): 



5t 

T 



5a 6 < V > 
-3.838 4.793 

a < V > 



(5) 



where we have used = —2^^. For the thermonuclear react ion rates dependence on a. 



we take the phenomenological expressions of tables IV and V in 



Landau et al. 



(120061 ). 



Since BBN is very sensitive to en, we should go further than equation (jlj) to evaluate 
the changes in in terms of 6 < v >. We will give our expressions for the different stages 
in terms of and ^y^. In section 4 we find an estimate for ^ in terms of 



<V> ' ££) ' O 



<v> 



and thus we can give then the final expressions in terms of 



5<v> 5a 
<v> ' a 



and 



5Q.Bh? 
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3. Abundances and their dependence on a and < v > in the different stages 

The ratio X of the number of neutrons to the total number of baryons in the first 
stage until the freeze-o ut of weak interactions ( T > 9.1 x lO^K) can be expressed as follows 



(IBernstein et al 



19891): 



1 / A2 







X{y = oo)= I dyey^-^{y) e^^HV = 0.151 (6) 



where K{y) = 6 + ^ + i + + ^) e"^] . Only b = 255^ depends on the 

fundamental constants through r and Am, so we get: 

X{y = oo) b 
Using equations [3] and [5] we obtain: 

^4P^^ = 1.385^-0.842^^ (8) 

X{y = oo) a <v> 

In the second stage, after weak interactions freeze out, neutrons decay freely until the 
rate of production of "^He becomes efficient (9.1 x lO^K > T > 0.93 x lO^K). Thus we have: 

y„ = X{y = oo) e"*/" = 0.151e-°-2/^9 (9) 

where Tg is the temperature evaluated in units of lO^K. 

The abundance of deuterium follows its equilibrium value and we assume the reactions 
[npd'-f] and [d'ynp] dominate for its production and destruction. Taking Qsh"^ = 0.0223 we 
can calculate the final temperature of this stage by setting 1^ = and thus 2Yn = — ^nl'^]- 
We find Tg = 0.93 and get the abundances Yp = 0.76 and F„ = 0.12. In order to calculate 
the dependence of the final temperature on the fundamental constants, we derive the 
equation 2y„ = — y„[n] with respect to the fundamental constants to get: 

— ^ = 0.065 ^ , „ + 0.055 0.119 + 1.195— 

Ig UbIi^ OL < V > tD 
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[npd'y] V , , , „, 1 



2^"(^P x^ rHi ) [ddpT] = - (10) 



and thus we get the relative variations of the nucleons abundances for this stage as follows: 

— = 0.030—4^ + 2.159 2.005 + 0.553— 

Yn ilBh^ a < V > ec 

— = -0.009^-0.682^ + 0.634^^-0.174^ 

In order to compute the final abundance of helium, we notice that once ^He production 
becomes efficient (i.e. 21^ = y„[n]), neutrons combine to form a-particles, and the 
production of the latter is dominated by [dTna] and [pT'^a]. One gets for the temperature 
of the ^ife freeze-out the following equation: 

2 

where Yp = 0.76, y„ = 0.151e-o-2/^9 and r is the neutron decay constant. Numerically 
we find Tg = 0.915 which is lower than the final temperature of the previous stage and 
larger than the final temperature of the next one. For the final helium abundance we find 
y/ = 2Yn = 0.238. As before, deriving equation [TO] with respect to e^, < > and a we find 

= 0.061— 4ir + 0.049 0.113 + 1.149— 

Tg" VLbH^ a <v> tD 

Since Yl = 2Yn we get the relative variation of the helium abundance as: 

5Yj dVlBh^ 5a 5 < v > 5en 
-J^ = 0.029— ^ + 2.182 2.042 + 0.549— 

YA ^sh? a <v> to 

In the following 'neutron cooking' stage, corresponding to 0.93 x lO^K > T > 
0.766 X lO^K, the neutron abundance can be expressed as follows: 

1 „ /■* f,. [npd-i] y ^ ~^ 



where the initial condition is given by the final values of the previous stage: Y^ = 0.12 at 
TgO = 0.93. 
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Putting Yn = Yd as the condition which determines the final temperature of this stage, 
we find y"^^^p] = 1- With Yp freezed at 0.76, we get numerically: 

T/ = 0.766 (12) 
F„ = 6.4xlO-^ = Frf (13) 

Again, the condition {Yn = Y^) allows the calculation of the relative change of the final 
temperature: 

— ^ = 0.031 ^ - 0.021— + 0.020 + 1.041 — 

and we get numerically: 

5Yd SYn SVLnh^ 5a 5 <v > Sen 
— = — = -1.095—4^ + 1.865 0.075 2.275— 

Yd Yn iiBh^ Oi <v> en 

In the last stage (T < 0.766 x lO^K), we notice that the dominant term in the time 
derivative of Yd is the production of Tritium, i.e. YdYd[ddpT]. We find: 

= (70 + 2 /* [ddpT]dt\ (14) 

\ d ^ ^initial ' 

with the initial value Y^ = 6.4 x 10"^ at = 0.766. 

We can obtain the final abundance of deuterium by setting the temperature equal to 
zero and we find the abundance numerically equal to yJ = 2.41 x 10"^. Again, we can 
numerically evaluate the relative change of yJ: 

^ = -1.072^ + 2.318^-0.049^^- 2.469^ 

YJ VLbH^ a <v> en 

We summarize the results that we obtained in table O 



4. The dependence of deuterium binding energy on the Higgs vev < v > 



As we said before, the deuterium binding energy e/5 is the most significant 
factor that can infiuence the BBN reactions rates, and its variation was discussed in 
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Table 2: Final abundances and their relative variations in terms of the relative variations of 



the fundamental constants, 



y/ 



<V> 



Y/ 


A 


B 


C 


D 


Abundance 


D 


-1.072 


2.318 


-0.049 


-2.469 


2.41 X 10-5 




0.029 


2.182 


-2.042 


0.549 


0.238 



Flambaum and ShuryakI (12002 



20031); 



Dmitriev and FlambaumI ((20031); 



Dmitriev et al. 



(I2OO4J ). Indeed, the equilibrium concentration of deuterons and the inverse reaction 
rate depend exponentially on it. Moreover, the deuterium is a shallow bound level 
(ez) ~ 2.225MeV). Therefore the relative variation of the deuterium binding energy en 
is much larger than the relative variation of the strong interaction potential which we 
neglect in our work. In order to give an estimate for the relative variation of eo, we should 
compute, within perturbation theory, the correction to €£> due to the perturbation which 
might change in time. Thus we write = ^% + AE, where e^) is the unperturbed binding 
energy and we co nsider it a QCD-deter mined quantity which does not change in time. As 



to AE, we know (jWeinberg 



1990 



1991 



that the one-pion exchange potential represents the 
first approximation to the perturbation on the strong interaction potential, and it has the 
form: 



V 



Y 



P e- 



An r 



(15) 



r2 

where ^ ~ 0.08, ~ 140MeV is the pion mass. We simplify the strong interaction 
potential by a square well model with width a and depth Vq. These two parameters 
can be determined by fitting the square well 'theoretical' expressions involving these 
two paramete rs to the p-ii scattering data. According to the shape-independent effective 



range theory (IBethe 



1949 



Bethe and Longmirei ll950l ) all the binding and low energy 
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scattering properties of the potential are determined by just two parameters which can be 
determined experimentally: the scatter ing length at = 5.50 x lO^^'^cm and the effective 



range Vt = 1.72 x 10 ^^cm (jSchifl 



19681 ). The corresponding values for the square well are: 
the depth Vq = 35.5MeV and the width a = 2.03 x lO^^^cm = 0.0103MeV"^ We consider 
the width and the depth as QCD-determined parameters and assume they do not change in 
time. 

Now, we have in the square well model 

AE = ——r / e e dx 



^J^l -^—"^ (16) 



where: 



^ = V^NiVo-eD) = 176.76 MeV (17) 



p = ^/^^^^ = 45.71 MeV (18) 

^ \ a sin{2^a) . sin'^{^a) 1393 McV (19) 

V 2 4^ ' 2/3 

B = -|sin(^a)e^"A = -559 MeV^/^ (20) 

and is the reduced mass of the two-nucleon system. 
Whence, 

AE = 0.203 MeV (21) 

On the other hand, the change due to the variation of the nucleon mass in the 
Yukawa potential V"^ is negligible compared with the change due to the variation of m^r. 
Since = as mentioned before, we get, evaluating numerically the integrals, the 

following: 

SAE 5mT, 5 < V > , , 

= -0.896 — - = -0.448 22 

AE < V > 
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and so, 



5e 



D 



-0.041- 



5 <v > 



CD <V > 

'hese values are one order of magnitude lower than those obtained by 



(23) 



Flambaum and Shuryak 



(|2002| ). The difference arises from the fact that these authors did not consider the continuity 
of the wave function and its derivative on the boundary of the square well. This results in 
differences between the normalization factor of the wave functions which propagate into the 
binding energy first order perturbation. Moreover, the values of a and Vq they consider are 
different from this work. 



Hence, the final expressions for the relative variations of the helium and deuterium 
abundances are: 



Yi 



-1.072—4^ + 2.318— + 0.052 



0.029 



5nBh^ 



6a 6 < V > 
2.182 2.044 

a < V > 



These results are summarized in table [3] 



(24) 
(25) 



Table 3: Abundances and their dependence on fundamental constants, f 

i 



f 



5<v> 
<v> 



Y/ 


Abundance 


A 


B 


C 




2.41 X 10-5 


2.318 


0.052 


-1.072 




0.238 


2.182 


-2.044 


0.029 
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5. Results 

We can now compare the theoretical predictions of the abundances of ^He and D 
obtained in the last section with the observational data. The equations fl24t [251) are of the 
form {i = l{D),2{^He)) 



5Yl_ 
Yl 



A, 



— — + B — + d 



sy/ 



(26) 



and we take the assumption that the difference is due to a change in the considered 



fundamental constants: ^ 



Yobs_YSBBN c R R h 

— ' ySBBN — , where and Y° are the theoretical and 



observed abundances respectively. In table HI the theoretical abundances Yf^^'^ are given 
for neh^ = 0.0223 with their errors resulting from the uncertainty in the values of the 
parameters involved. In table [5], the observational data of helium and deuterium are stated 
with their measured errors. 



Table 4: Theoretical abundances in the standard model with the WMAP estimate ils/i^ 
0.0223 



Nucleus 


ySBBN _|_ ^ySBBAf 




(2.51 ±0.37) X 10^5 


^He 


0.2483 ±0.0021 



As regards the 



Landau et al 



consistency of the D and "^He data, we follow the treatment of 



([20061) and increase the observational error by a factor O. The values of G 



are 2.4 for D, 2.33 for ^He. 



The results of solving the system of equations f[26|) with the given data are shown in 
table [6l These results are consistent within 1 — a with no variation of the fundamental 
constants. On the other hand, the results considering variation of one fundamental constant 
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Table 5: Observational abundances 



Nucleus 


yobs ^ ^Y°^s 


Reference 


D 


(1.65 ±0.35) X 10-5 


Pettini and Bowen f2001) 


D 


(2.54 ±0.23) X 10-5 


O'Meara et al. f2001) 


D 


(2.42l°;i) X 10-5 


Kirkman et al. f2003) 


D 


(3.25 ±0.3) X 10-5 


Buries and Tvtler fl998al 


D 


(3.98t[J:59) X 10-5 


Buries and Tvtler f 1998b) 


D 


(1.61^:^5) ^ 10-5 


Criehton et al. f2004) 


^He 


0.244 ±0.002 


Izotov and Thuan fl998) 


^He 


0.243 ± 0.003 


Izotov et al. ri997l 


^He 


0.2345 ±0.0026 


Peimbert et al. f2000) 


^He 


0.232 ±0.003 


Olive and Steieman f 19951 



only are shown in table [3 These results are consistent within 3 — cr with no variation of the 
fundamental constants. 

Fig. 1. — The full line shows the theoretical probability of the residuals, while the dotted 
line shows the empirical probability. Deviations of Vtsh'^^ oi and < v > with respect to their 
mean values are considered 

0.9 
0.8 
0-7 
0.6 
Q- 0.5 
0.4 
0.3 
0.2 
0.1 
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Table 6: Results of the analysis using the square well model considering joint variation of all 
constants 



Realtive variation 


Value 


a 


SQBh^ 


-0.011 


0.054 


5a 
a 


-0.032 


0.072 


5<v> 
<v> 


-0.011 


0.078 



Table 7: Results of the analysis using the square well model considering variation of each 
constant only 



Realtive variation 


Value 


a 


SUBh^ 

nBh? 


-0.013 


0.056 


&a 
a 


-0.015 


0.006 


S<v> 
<v> 


0.017 


0.007 



In order to check the goodness of our fit, we performed a Kolmogorov-Smirnov (K-S) 
test (see figure E]). When considering variation in fi^/i^, a and < v > altogether, we have a 
probability of 79% to obtain a worse fit. However, we consider the results of the K-S test 
only indicative, since even though the considered data are independent the residuals are 
not. 



Finally, it is interesting to compare our results with other non-standard BBN models. 
In particular, a non-standard expansion rate and an electron-neutrino a symmetry were 



considered by several authors (IBarger et al. 



2003a 



Steigman 



2005 



20061 ). While D is more 



sensitive to the baryon density [Qfyh'^), the effect of a non zero z/e asymmetry is more strong 
for ^He than for the other relic nuclides. This is similar to the effect of < f > as follows 
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from equations [23] and [25], whereas both D and ^He are sensible to changes in a. In the 
papers cited above, the authors considered the following cases: i) adding only one free 
parameter (either AA'^^ or ^g) and ii) adding two free parameters [AN^, and C,e) besides 
the baryon density. They consider the WMAP data from the CMBR togheter with ^He 
and D abundances. In the first case, the results are consistent within 1 — a with 7^ 
and AN^, 7^ 0. In the second case, the results are consistent within 1 — a with = and 
AA^j, = 0. These results are similar to ours in that considering only one free parameter 
besides the baryon density is consistent with non standard physics within 1 — a whereas 
considering two free parameters is consistent within 1 — a with the SBBN model. However, 
it is important to remind that most of the theories where the fundamental constants may 
vary in cosmological time scales, predict joint variation of constants. 

6. Conclusions 

In this work, we assumed that the discrepancy between SBBN estimation for '^He and 
D and their observational data is due to a change in time for the fundamental constants: 
the Higgs vev < v >, the fine structure constant a. We analyzed the dependence of the 
'^He and D abundances on these fundamental constants within perturbation theory and 
on deviations with respect to the mean value of the baryonic density. Furthermore, we 
compared them with the observational data. We find that varying fundamental constants 
may not solve, in our case, the discrepancy between the theoretical SBBN and the observed 
data considered in this work. We hope this work stimulates further research in this 
interesting subject. 
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